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TRANSIENT  AND  LOAD-FOLLOWING  CHARACTERISTICS  OF  A  FULLY 
INTEGRATED  SINGLE-CELL  THERMIONIC  FUEL  ELEMENT 

ABSTRACT 

Practical  Application  of  thermionic  space  nuclear  power  systems  necessitates  the 
investigation  of  such  systems  le^nses  to  changes  in  operation  and  design  parameters  as  well  as 
load-following  characteristics.  The  principle  huilding  block  in  these  systems  is  the  thermionic  fuel 
element  (TFE).  A  Thermionic  Transient  Analysis  Model  (TITAM)  is  developed  to  simulates 
tranaent  and  steady-state  operations  of  a  fiilly  integrated,  single-cell  TFE.  TITAM  is  used  to 
investigate  the  responses  of  the  TFE  to  a  step  input  in  reactivity,  a  change  in  Cs  pressure  and/or  in 
the  size  of  the  interelectrode  gap,  a  change  in  the  coolant  temperature,  and  a  change  in  load 
demand.  The  effects  of  these  parameters  on  load  electric  power,  emitter  temperature,  overall 
conversion  efficiency,  and  load-following  characteristics  of  die  TFE  are  determined.  Results  show 
that  although  nuclear  reactors  having  negative  temperature  reactivity  coefficients  are  always  load¬ 
following,  TFEs  are  only  partially  load-following.  However,  for  a  given  load  electric  power  need 
there  are  several  combinations  of  fission  powers  and  load  residences.  Results  also  show  that  for 
TFEs  having  a  large  interelectrode  gap,  it  is  desirable  to  conserve  Cs  by  lowering  its  vapor 
pressure  at  the  beginning-of-life,  since  increasing  the  Cs  pressure  insignificantly  affects  the  load 
electric  power.  However,  should  fuel  swelling,  after  operating  the  reactor  for  an  extended  period 
of  time,  reduces  the  width  of  the  interelectrode  gap,  both  the  cemversion  efficiency  and  the  load 
electric  power  will  decrease.  In  this  case,  the  load  electric  power  could  be  restored  by  increa^g 
die  fission  power,  and  only  partially  by  changing  the  coolant  temperature  and/or  increasing  the  Cs 


vapor  pressure. 
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Ouqnt  voltage  at  the  transitioo  point  in 
the  ignited  mode  (V) 

Emitter  Aeatb  voltage  drop  (V) 
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a 
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ft 
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ft 
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5 
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♦ 
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c 
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f 
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m 
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n 
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1.  INTRODUCTION 


Several  space  exploration  missions  are  being  considered  early  next  century  which  would 
require  electric  power  levels  ranging  from  tens  to  hundred  of  kilowatts.  Examples  of  such 
missions  ffiff^hiishing  an  outpost  on  the  moon  and  mining  of  its  resources,^  ^  powering  a 

manned  rovers  for  global  e]q>loration  of  the  Lunar  and  Martian  surfaces, and  powering  a 
mynnfiH  locket  to  Mars.®"®  The  power  needed  for  these  and  other  defense  and  civilian  missions 
operating  in  high  earth  orbit  for  3-10  years  are  best  met  with  the  use  of  nuclear  reactor  ^ace  power 
systems.  Deagn  requirements  of  such  systems  include  safety  during  all  phases  of  the  mission 
(launch,  deployment  and  end-of-life  disposal),  low  specific  mass,  no  single  point  failure, 
survivability,  low  heat  rejection  radiator  area,  and  hi^  reliability. 

Nuclear  reactor  power  systems  are  currently  being  developed  under  two  joint 
DOE/DOD/NASA  programs  to  advance  the  Thermionic  (TI)  and  Thermoelectric  (TE)  power 
systems  technologies^*  to  support  a  potential  nuclear  flight  demonstration  program  decision  in 
mid  1990s.  Both  systems  offer  many  attributes  for  meeting  electric  power  needs  from  a  few  tens 
to  hundreds  of  kilowatts.  The  TI  nuclear  reactor  power  systems  have  a  moderately  high 
conversion  efficiency  (8%-10%)  and  operate  at  higher  radiator  temperature  (>900  K)  than  TE 
systems.  The  high  efficiency  lowers  the  reactor  thermal  power,  resulting  in  lower  masses  for  the 
reactor  and  shield  and  smaller  radiator  area.  Also,  the  high  radiator  temperature  reduces  its  size, 
resulting  in  smaller  system’s  overall  volume  and  specific  mass.  Like  its  sister  TE  systems,  TI 
systems  promote  high  reliability,  due  to  the  lack  of  moving  parts,  and  redundancy,  since  a  failure 
of  a  single  converter  unit  will  have  little  effect  rai  Ae  performance  of  odier  unites  in  the  system. 
However,  since  TI  converters  could  be  used  boA  in-core  or  ex<ore,  Aey  provide  the  system 
designer  a  flexibility  in  integrating  Ae  system  components  for  oAanced  overall  reliability  at  a 
potratial  reduction  in  tiie  overall  system  mass  and  volume. 

An  important  issue  related  to  Ae  development  of  TI  q>ace  nuclear  power  systems  is 
understanding  Ae  effects  of  various  operation  and  design  parameters  on  Ae  steady-state  and 
transient  performance  of  Ae  fully  mtegrated  power  system,  or  Ae  Aermionic  fuel  elements  (TFEs), 
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which  are  the  basic  building  blocks.  Some  of  these  parameters  include  changes  in  reactor  thermal 
power  to  a  reactivity  insertion,  changes  in  the  width  of  interelectrode  gap  due  to  fuel  swelling, 
changes  in  Cesium  (Cs)  vapor  pressure  in  the  interelectrode  gap,  and  changes  in  the  load  demand, 
b  order  to  quantify  the  effects  of  these  parameters,  an  integrated  system  or  TFE  model  needs  to  be 
developed,  which  directly  coiq>les  die  TI  converter  umts  to  both  the  electric  load  and  die  nuclear 
fuel  region  and  incorporates  nuclear  reactor  kinetics  widi  the  appropriate  reactivity  feedbacks. 

Existing  models  are  limited  to  die  analysis  of  steady-state  operation  of  either  a  TI  diode**  or 
aTFE.*2-*3  The  input  to  the  TI  diode  model  includes  die  emitter  and  collector  temperatures  and  the 
emission  current  density.  The  TFE  model  of  Pawlowski  et  al.*2  basicaUy  couples  the  TI  diode 
model  to  a  two-dimentional,  steady-state  thermal-hydraulic  model.  Because  this  model  is  not 
directly  coupled  to  the  electric  load  and  reactor  kinetics,  the  load  current  and  fission  power  must  be 
specified  in  the  input.  In  addition,  both  models****^  do  not  include  the  unignited  mode  of  TI 
emission,  hence,  are  not  suited  for  startup  simulation. 

The  objective  of  the  research  is  to  develop  a  transient  model  of  a  fully  integrated  single-cell 
TFE,  which  directly  couples  the  electric  load  to  the  TI  converter  and  incorporates  neutron  kinetics 
with  the  appropriate  reactivity  feedback  effects.  The  model  is  used  to  investigate  the  effects  of  a 
reactivity  step-function  input,  a  change  of  the  Cs  vapor  pressure  and/or  the  width  of  the 
interelectrode  gap,  and  a  change  in  the  coolant  temperature  <®  tiie  TI  conversion  efficiency,  load 
electric  power  and  TI  converter  temperatures.  The  load-following  characteristics  are  also 
investigated  and  tiie  operation  regimes  in  which  tite  TFE  is  inherently  load-following  are  identified 
as  functions  of  fission  power  and  load  demand. 

2.  MODEL  DESCRIPTION 

The  Thermionic  Transient  Analysis  Model  (TITAM)  is  developed  to  simulate  transient  and 
steady-state  operations  of  a  fully  integrated  TFE  uith  eitiier  fission  heating  or  electric  heating.  The 
fission  heating  option  could  be  used  to  predict  the  performance  chararteristics  of  in-core  TFEs, 
while  the  electric  heating  option  would  be  useful  in  comparing  the  model  predictions  with  test 
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results  of  electrically  heated  TFEs.  THAM  consists  of  four  interactively  coupled  sub-models:  (1)  a 
nuclear  reactor  six  groups  point-kinetics  model,  (2)  a  lumped  transient  thermal  model  of  the 
integrated  nuclear  fuel  element  and  TI  converter,  (3)  a  thermal-hydraulics  model  of  the  liquid  metal 
coolant  and  (4)  a  thermionic  emission  model.  **  The  woric  functions  of  the  emitter  and  the 

collector  materials  in  die  TI  emission  sub-model  are  determined  as  functions  of  the  electrode-Cs 
reservoir  temperature  ratio  u^g  cwrelations  based  on  the  low  Cs  pressure  data  of  Houston, 
Ongoing  effort  is  directed  toward  updating  tiiese  correlations  based  on  hi^  pressure  data  (up  to  20 
torr)  and/or  jdienomenological  modeling. 

2. 1  Input  and  Output  Parameters 

In  die  fission  heating  option,  the  input  to  TITAM  includes  the  dimensions  and  materials  of 
die  TFE,  temperature  reactivity  feedback  coefficients  of  the  different  regions  of  the  TFE,  cesium 
reservoir  temperature,  load  electric  resistance,  either  fission  power  or  reactivity  input  to  the  reactor 
core,  and  the  coolant  inlet  temperature  and  mass  flow  rate.  The  output  parameters  includes  electric 
power  to  the  load,  load  current  and  voltage,  conversion  efficiency,  and  temperatures  in  the 
different  regions  of  the  TFE.  For  electrically  heated  TFE,  the  input  and  output  parameters  are  the 
same  as  for  the  fisaon  heating  option,  except  that  fission  power/reactivity  input  is  replaced  by 
electric  power  input,  with  no  reactivity  feedback  effects.  As  a  user  option  the  dimenaons  of  the 
electric  leads  could  be  optimized  for  minimum  electric  and  thermal  losses.  Figure  1  presents  a 
and  the  electric  circuit  of  the  single-cell  TFE  model  and  Fig.  2  presents  a  line  diagram  of 
die  building  blocks  in  TITAM. 

2.2  Capabilities  of  TITAM 

hi  addition  to  transient  and  steady-state  simulation,  TITAM  could  be  used  for  design  and 
operation  optimization  of  single  cell  TFEs  and/or  of  a  TI  power  system.  An  example  of  potential 
in-core  application  of  TITAM  is  mapping  the  design  parameters  in  the  nuclear  reactor  core  for 
maintaining  a  coiistant  emitter  temperature  and  assessing  die  effects  of  non-uniform  radial  and  axial 
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Figure  1.  A  Schematic  of  the  Single-cell  Thermionic  Fuel  Elemrat 
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Figure  2.  A  Line  Diagram  of  die  Building  Blocks  in  TII'AM. 


distributions  of  rission  power  in  the  reactor  core  on  the  performance  of  the  TFEs.  Examples  of 
other  {plications  include:  (a)  determining  the  operation  regime  udiere  the  TFEs  or  the  TI  power 
system  is  inherently  load-following,  (b)  determining  the  effect  of  using  high  temperature  heat  pipe 
to  equalize  the  temperature  almg  the  emitter  on  die  TFE  performance,  (c)  developing  an  operation 
strategy  for  ptimizing  the  cmisumption  of  Cs  over  the  lifetime  of  the  system,  and  (d)  determining 
the  effects  of  using  a  low  pressure  Cs-Ba  vapor  mixture  (  <  0. 1  torr),  instead  of  die  hig^  pressure 
Cs  vapor  (2-10  torr),  on  the  performance  of  a  new  generation  of  high  performance  TFEs. 

The  TI  sub-model  in  TTTAM  incorporates  both  the  ignited  and  the  unignited  modes  of 
emission,  uriiich  makes  it  suitable  for  both  high  temperature  and  low  temperature  simulations 
including  startup  and  shutdown  transients.  Because  the  work  functions  of  the  electrode  are 
ttetermined  in  terms  of  the  Cs  reservoir  temperature  and  those  of  the  electrodes,  TIT AM  could  be 
used  to  investigate  the  effects  of  changing  the  Cs  pressure  on  the  load  electric  power,  overall 
converrion  efficiency,  and  the  load-following  characteristics  of  the  TFE.  As  indicated  earlier,  in 
addition  to  simulating  the  operation  of  a  fission-heated  angle  cell  TFE,  by  including  temperature 
and  Doppler  reactivity  feedback  effects,  TIT AM  is  capable  of  simulating  electrically  heated  TFEs 
uhich  could  be  tested  in  a  laboratory  setup.  Therefore,  the  user  would  be  able  to  relate  the 
performance  of  electrically-heated  TFEs  to  that  of  fission  heated  TFEs,  without  the  expense  and 
complicaricm  of  fission-heated  tests.  The  component  models  shown  in  Hg.  2  ate  discussed  in  detail 
in  the  following  secticxis. 


2.3  Six  Groups  Point-Kinetics  Model 

The  six  groups  point-kinetics  model  in  TIT  AM  calculates  the  transient  fission  thermal 
power  following  an  external  reactivity  insertion  and  incorporates  reactivity  feedback  doe  to  Doppler 
and  temperature  effects.  The  nuclear  reactor  point-kinetics  equations  may  be  expressed  as : 


dPinit)_pit) 

’  A 


^»i(0  + 


£  A/C;(r) 
r^l  ,and. 


(1) 
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The  total  reactivity  input  to  the  reactor  core  is  g^ven  as: 


P(0  =  P«(0+p/(0  +  2pf(0 

The  first  term  on  flie  R.RS  of  Eq.  (3)  is  the  external  reactivity  input  to  the  nuclear  reaaor  core. 
The  second  and  third  terms  in  Eq.  (3)  are  the  reactivity  feedback  due  to  Doppler  effect  in  the  fuel 
region,  PpW,  andPf^^)  is  die  gim  of  temperature  reactivity  feedback  effects  in  different  regions  of 
the  TFE.  The  Doppler  feedback  reactivity  and  the  temperature  feedback  reactivity  in  any  given 
region,  i,  in  die  TFE  [see  Fig.  1]  are  expressed,  respectively,  as: 

{Jf{t)rrf{0)\  (4) 

and, 

pf(r)  =  af(rK()-77«»).  (5) 

2.4  Thermal  Model  of  the  TFE 

The  thermal  sub-model  in  TITAM  calculates  the  average  transient  and  steady-state 
temperatures  in  the  various  regions  of  the  TFE  [see  Fig.  1]  including  the  emitter,  collector, 
insulation,  cladding,  and  liquid-metal  coolant  in  the  cooling  channel,  as  functions  of  fission  power 
(or  reactivity  input),  coolant  flow  rate  and  iidet  temperature,  width  of  interelectrode  spacing,  Cs 
vapor  pressure,  and  other  design  parameters.  The  average  fuel  and  fuel/emitter  gap  temperatures 
during  a  tranaent  are  given,  respectively,  as: 


—  Pin  •  hfAfiTf^  Tgaf 

MgcPga^  =  hfAfiJf-  Tg^  -  hgaAga{Tga  -  Tjf^ 


where  hf^  (2  kfirfh  For  an  open  fuel/emitter  gap  the  gap  conductance  is  calculated  as: 


3aTh 


Sgal  + 


where  =  ^(Tfe)  Tfe  =  VTfTe.  The  average  emitter  temperature  and  that  of  the  interelectrode 
gap  are  given,  respectively,  as: 

^  S  hgaAga  (Tga  -  Te)  -  M  e(Tga  -  Te)  •  hisA  UTe  -  Tc)  -  Pec  +  (9) 


MgCg—^^  heA^Te-Tg)-  hgAg{Tg-Tc) 

In  Eq.  (10),  the  radiation  and  conduction  contributions  to  heat  transfer  between  the  emitter  and  the 
collector  are  incorporated  in  the  gap  conductance,  hg,  as: 


'  (I)  *  (ife  +  ite-i 


In  this  equation,  -f  (Te  )♦  Be  =/  (Tec  )♦  Tgc  =  The  thermal  conductivity  of 

Cs  vapor  in  the  interelectrode  gap,  kg,  is  determined  in  terms  of  its  vjqwr  pressure  as:** 


1.65x10-^<TZ 
1,14.  10-5 _ 


(12a) 


where. 
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The  tran^ot  temperatures  of  tiie  collector,  insulator  and  cladding  are  givra,  respectively,  as: 


hgAf^Tg  kcA(^Tc -  Tins)  ^  hitAi^Tg~Tc)  +  Pec ^2) 

®  f^cAc  iXc  "  ^ins)  "  f^insAinsiPins  "  Td) 

and, 

MdCcl  *  hins  Ain^Tins  -  Td)  *  hdAclfJd  *  Pco)  ( j5) 

2.5  Coolant  Channel  Thermal-Hydraulics  Model 

The  transient  bulk  temperature  of  the  liquid-metal  coolant  in  the  annular  flow  channel  of  die 
TFE  [see  Rg.  1]  is  determined  from  the  relation: 

^coAalTd  ~Tco)-  2tHCco{Tco~Tin) 

uhere  the  heat  transfer  coefficient  for  concentric  annuli  is  given  as:^^ 

hco^{kcoff>e)iA*BPe^).  (17) 

where, 

i4  *4.63 +  0.686  y,  (18a) 

B  =  0.215  +  4.3x10-5  y,  (18b) 

c  *0.752  +  0.01657  y- 8.83x10-4)2.  (18c) 


During  steady-state  simulation  the  terms  on  the  left-hand  side  of  Eqs.  (6),  (7),  (9),  (10), 
(13),  (14),  (15)  and  (16)  are  dropped  out  and  the  temperatures  in  the  different  regions  of  the  TFE 
are  calculated  using  an  iterative  explicit  approach,  hi  fliis  case,  the  iterations  are  only  for  the  sole 
purpose  of  updating  tiie  thermophyacal  properties  in  Ae  various  reports  as  functions  of  their 
respective  temperatures. 

2.6  Thermionic  Emission  Model 

The  structure  of  the  TI  emission  model  in  TUAM  follows  closely  the  phenomenology 
developed  by  Rasor,i^  ^hich  takes  into  account  the  continuity  of  electron  and  ion  emission 
currents,  and  energy  fluxes  at  the  plasma-electrode  interfaces  and  across  the  plasma.  The  complex 
collisional  effects  in  the  plasma  have  been  neglected.  The  thermionic  emission  model  calculates  the 
electric  power  delivered  to  the  load  in  a  self-consistent  manner  with  the  thermal  sub-model 
discussed  earlier.  The  TI  emission  model  include  both  the  ignited  and  unignited  modes  of 
operation.  Also,  the  transition  conditions  between  the  saturated  and  obstructed  regimes  in  each 
mode  are  included  in  the  models.  The  overall  energy  balance  in  the  thermionic  emission  model 
includes  electron  cooling,  heat  conduction  through  the  interelectrode  gap,  and  heat  conduction  and 
Joule  heating  in  die  electric  leads.  The  characteristic  equations  describing  these  processes  are: 

Electron  Cooling 

+  (19) 

Heat  Conduction  in  the  interelcctrode  gap 


(20) 


Heat  Conduction  through  the  leads 


^de  -  (^e  “  ^c) . 


(21) 


Joule  heating  in  the  leads 


Pjle^iJAefRie,  (22) 

A  steady-state  the  overall  eneigy  balance  for  the  TI  converter  yields : 


out  ^  Pload  ,  where, 

(23) 

Pin^  Pcle"  0.5  Pjle^ 

(24a) 

Pload^  ^t^^loady 

(24b) 

and, 

^ load  out  die^ 

(25) 

A  quasi-steady  state  approach  is  used  to  couple  Eq.  (23)  with  the  transient  thermal  model 
alien  simulating  the  transient  operation  of  the  TI  converter.  During  each  small  time  step  in  die 
diermal  model  it  is  assumed  that  the  condititHis  governing  thermionic  emisaon  are  constant,  but 
different  from  those  in  die  previous  and  subsequent  time  steps.  This  approach  should  be  fairly 
accurate  when  die  time  step  is  small  raou^,  depending  of  the  sdffiiess  of  the  problem,  hi  fact, 
because  the  sdfhiess  of  die  reactor  kinetics  equations,  the  time  step  used  for  tianaent  fission  power 
c^iniiatinns  ( ~  0. 1  sec),  easily  satisfied  the  requirements  for  die  quasi-steady  state  soluticm  of  the 
TI  emission  model.  This  conclusion  is  reached  vdien  it  is  found  that  a  decrease  in  die  time  step. 


beyond  dictated  by  tiie  leactor  kinetics  equations,  does  not  practically  change  the  results. 
Also,  to  ensure  numerical  stal^ity  of  Ae  transient  solution  and  flexibility  in  Ae  selection  of  Ae 
fimft  step,  an  implicit  time>Ategration  approach  is  used  to  solve  Ae  governing  equations  A  Ae 
Aermal  model  of  the  TFE. 


3.  THERMIONIC  EMISSION  MODEL 

TTie  emitter  and  collector  tempmtures,  cerium  pressure,  and  electrode  work  functions 
determine  Ae  operating  mode  of  Ae  TI  converter.  At  lower  temperatures,  Ae  converter  is  m  Ae 
unignited  mode  and  at  higher  temperatures  (>1300  K)  Ae  converter  is  m  Ae  ignited  mode.  In  each 
mode,  Ae  output  voltage  (V,„,)  is  determined  from  Ae  J-V  characteristics  for  Aat  mode.  The  TI 
emission  models  for  Ae  ignited  and  unignited  modes  of  operation  are  presented  m  Ae  following 
sub-sections. 

3.1  Ignited  Mode  of  Operation 

In  Ae  ignited  mode  of  operation,  Ae  melastic  collisions  m  the  mterelectrode  space  is  Ae 

<)ofninanf  ionization  mechanism  of  Cs  atoms  m  Ae  gap.  The  current  discharge  m  the  g^  is 

sustained  by  Ae  potential  Aop  across  Ae  mterelectrode  space,  Vj.  When  Cs  ions  generated  m  Ae 

gap  by  inelastic  collisions  are  just  sufficient  to  neutralize  Ae  ^ace  charge,  Ae  converter  is  said  to 

1  • 

be  operating  at  the  transition  pomt,  ’^out\  when: 

(26) 

The  transition  pomt  Aus  corresponds  to  Ae  mAimum  voltage  drop  required  to  fully 
neutralize  Ae  space  charge  for  a  given  P^gd,  resuItAg  A  a  maximum  conversion  efficiency.  The 
ffgiiim  vtqmr  pressure  A  Ae  Aterelectrode  gap  can  be  expressed  A  terms  of  Ae  cesium  reservoir 


temperature  as:*^ 


(27) 


Id  tiie  ignited  mode  of  operatira  dieie  are  two  distinct  regions  i^liich  intersect  at  the  transition 
point;  the  obstnicted  ignited  emission  region  (Vont  >  saturation  ignited  emission 

region  (Voot  <  ^o«).  These  two  regions  are  characterized  by  distinctly  different  processes,  hi  the 

obstnicted  region,  because  is  insufficient  to  fully  neutralize  the  space  charge  a  negative  space 
diarge  barrier  DV  forms  at  die  emitter.  This  potential  barrier  reduces  die  electron  emission  current, 
dius  making  space  charge  neutralization  possible,  hi  the  obstructed  ignited  region,  the  discharge 
current  is  given  as:^^ 


and  the  cnitput  voltage  is  given  as: 


Vcur  =  ^e-<l>c-Vd+AV  ^  (29) 

where  Je,  R,  T^g,  (jf),  are  defined  by  Eqs.  (9),  (24-29),  and  (31)  in  Ref.  14. 

In  the  saturation  ignited  region,  Vj  is  greater  than  which  is  required  to  completely 
neutralize  the  space  charge.  The  excess  potential,  -AV  ^(V in  the  emitter  riieath  generates 
positive  ions,  which  move  toward  the  emitter  giving  rise  to  a  positive  ion  current  The  electric  field 
created  at  the  emitter  surface  lowers  the  emitter’s  effective  worit  function,  hence  increasing  the 
electron  emissitm  by  Schottky  effect.  The  output  current  in  the  saturated  ignited  region  can  be 
expressed  as: 


J 


where  J'  is  the  coirent  density  at  the  transition  point,  J,  and  J,  are  given  by  Eqs.  (1),  (2)  and  (35) 
in  Ref.  (14).  Eqs.  (27-30)  are  solved  iteratively  to  determine  the  transition  point  and  generate  the 
J-V  curves  Dibich  are  characteristic  of  the  design,  electric  load,  and  operating  conditions  of  the  TI 
converter. 


3.2  Unignited  Mode  of  Operation 

In  the  unignited  mode  of  operation,  surface  ionization  is  the  dominant  mechanism  for 
maintaining  neutral  plasma  in  die  interelectrode  gap  of  the  converter,  hi  this  mode  of  operation, 
volume  ionization  is  absent  and  the  plasma  potential  can  be  expressed  as: 

^„*0.5[/i(r^/ti)+V,-ifc/»r^  (31) 


where  /i  =  0.76  eV  and  V,  =  3.89  eV.  The  electron  emission  current  is  given  as: 

,32) 


This  electron  current  is  partially  neutralized  by  the  ion  current,  which  is  given  by  the  Saha- 
Langmuir  equation  as:^^ 


/,  =  e/i 


1  +  2cxp 


\YiiM 

I  ksTe  /J 


-1 


(33) 


In  diis  equation,  the  ions  arrival  rate,  p,  is  given  as: 


li^Dexpi-hlksTR) 


(34) 


In  Eq.  (32),  A  =  120  A/m^  and  in  Eq.  (34),  D  is  taken  equal  to  (lO^V  Vl*)  atoms^cm^sec^l  fa 
the  unignited  mode  of  operation,  there  are  three  regions  of  interest:  (a)  retarding  region  (Vob,  > 
Vout).  (b)  transition  point,  (J",  Vont),  and  (c)  plasma  saturation  region  (Vont  <  Vo^)*  At  the 
transition  point,  (J".  Vo«),  which  coire^nds  to  the  condition  of  zero  electric  field  at  the  collector 
(no  collector  sheath),  the  current  density  and  voltage  relations  are : 


fa  the  retarding  region  (Vout  >  Von,)  tbe  current  density  is  given  as: 

# 

J-Jn  [(35^/4A  )  +  exp  ((Vout  -  ^out  \ 


(35) 


(36) 


fa  the  plasma  saturation  region  (Vout  <  V^^),  the  current  density  is  limited  by  the  ions  reflected 
from  the  collector  dieath  and  is  given  by: 


/  =  2/"-  Jn  [(35^/4A  )+ exp  ((V^,  -  Vom 


-1 


(37) 


Equations  (31-37)  are  solved  iteratively  to  determine  the  J-V  characteristics  of  the  TI  converter  in 
the  unignited  mode  for  various  design,  electric  load,  and  operating  conditions. 

At  the  transiticMi  betwem  the  ignited  and  unignited  modes  of  operation,  the  proper  operation 
regime  for  the  TI  cwiverter  is  determined  by  the  model,  depending  on  vdiether  the  TFE  is  in  a 
startup  or  a  shutdown  {diase.  During  foe  startup  of  foe  TFE,  foe  TI  converter  operates  initially  in 
foe  imignitftri  mode,  where  foe  emistion  current  increases  and  foe  arc  voltage  drop  decreases  as  foe 
fistion  input  power  is  increased.  When  foe  voltage  drop  in  foe  converter  decreases  below  that 
corresponding  to  the  intersect  between  the  ignited  and  unignited  characteristic  curves,  there  are  two 
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possible  values  of  the  emission  current,  the  low  value  describes  the  converter  operation  in  the 

unignited  mode,  and  the  higher  value  cotregionds  to  the  operation  in  the  ignited  mode.  In  this  case 

the  model  takes  the  higher  emission  current  and  declares  the  TI  converter  operating  in  the  ignited 

mode.  Beyond  this  point,  further  increase  in  fisaon  power  results  in  higher  mnission  current  and 

lower  voltage  drop,  however,  the  operation  regime  (saturated  or  obstructed)  is  determined  by 

,  • 

comparing  the  values  of  J  and  with  those  corresponding  to  die  transition  point  ♦  (J 
During  shutdown,  or  a  decrease  in  the  fission  power  of  the  TFE,  a  similar,  but  opposite  logic  is 
used  in  die  model  to  detomine  the  transition  from  the  ignited  to  the  unignited  mode  of  operation. 

3.3  Electrodes  Work  Functions  Correlations 

The  effective  work  function  of  the  electrode  depends  on  the  degree  of  cesium  coverage  of  its 
surface,  nhich  depends  i^on  the  electrode  surface  temperature  and  the  Cs  vapor  pressure  in  the 
interelectrode  gap.  Using  a  semi-empirical  curve-fits  of  low  Cs  pressure  experimental  data,  **  the 
work  functions  are  correlated  in  terms  of  the  ratio  of  the  temperature  of  the  emitter  or  collector  to 
the  cesium  reservoir  temperature  as  follows: 

Tungsten  Emitter 

)«1.683+1.69l(r,/rjf).2.646(rtf/rj?)?+  1.321(1^/1/})^- 0.244(re/rj?)*+0.0152(re/r/})*.  (38) 

Molybdenum  Collector 

(Mo)=1.1173+2.0169  (rc/r/}>1.7643  (Tc /!/})?+ 0.4405  {Tc/TRf  ,  fQit(Tc  /Tr)<  2.65^  (39a) 


9c(Wo)»-7.29+5.2^rc/r«>o.53(rc/r*)F+o.07{rc/r*)?  +o.02(rc/rii)*  -7.iuio^  {T/Tiif , 

for  (Tc/Tk)^  2.65 


(39b) 


These  conelaticms  are  compared  in  Rg.  3  with  the  low  pressure  data  of  reference  (15).  As  this 
figure  indicates  the  correlafions  fits  the  data  to  within  ±  5%. 

4.  RESULTS  AND  DISCUSSIONS 

The  steady-state  and  transent  analyses  presented  herein  are  for  optunized  electric  leads  and 
for  the  TFE  base  case  design  and  operation  parameters  listed  in  Table  1.  The  following 
subsections  present  steady-state,  transient,  and  load-following  results  for  a  fully  integrated, 
single-cell  TFE.  The  steady-state  results  show  the  effects  of  changing  the  fission  power,  Cs 
pressure  in  the  interelectrode  gap,  the  size  of  the  interelectrode  gap,  and  the  coolant  temperature  on 
the  TFE’s  load  electric  power,  temperatures,  and  conversiOT  efficiency.  The  transient  results 
the  response  of  the  TFE  to  both  a  step-function  increase  and  a  step-function  decrease  in 
reactivity.  Transient  results  include  the  fission  power,  load  electric  power,  and  the  TFE 
temperatures  a"d  the  overall  conversion  efficiency.  The  overall  converoon  efficiency  of  the  TFE  is 
defined  herein  as  the  ratio  of  the  load  electric  power  to  the  actual  fission  power 
generated  in  the  TFE  during  the  transient.  At  steady-state,  the  fission  power  is  equal  to  the 
riiermal  power  supplied  to  the  emitter.  However,  during  a  transent  increase  or  decrease  in  fission 
power,  the  corresponding  thermal  power  to  the  emitter  is,  respectively,  lower  or  higher  than  the 
fission  power.  Consequently,  the  TFE  conversion  efficiency  during  a  transient  would  be  markedly 
differem  than  the  TI  converter  efficiency.  The  load-following  results  define  the  operation  regimes 
and  the  ranges  of  operating  parameters  for  ahich  the  TFE  is  inherently  load-following  and  diose 
causing  the  TFE  to  be  ntmload-following. 

4.1  Steady-State  Operation 

Rgure  4  presents  die  effects  of  fisriem  power,  Cs  pressure,  and  aze  of  the  interelectrode  gtq> 
on  steady-state  electric  power  to  the  load  and  the  TFE  conversion  efficiency.  Results  show  that  for 
the  smaller  gap  size  (0.094  mm)  increasing  the  Cs  pressure  generaUy  results  in  both  higher  electric 
power  to  the  load  and  higher  conversion  efficiency,  however,  the  effect  is  more  pronounced  at 
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Hgure  3.  Woik  Functions  Conelations  for  Tungsten  and  Molybdenum  Electrodes  Based  on  the 
Data  of  Houston 
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Table  1.  Base  Case  Deagn  and  Operation  Parameters  of  the  TFE 


Parameter 


Value 


TICeUHei^t(cm) 

Fiiel  Gap  Size  (mm) 
Interelectrode  Gap  width  (mm) 
Emitter  Thickness  (cm) 


Collector  Thickness  (cm) 


Insulation  Thickness  (cm 


Cladding  Thickness  (cm) 
Coolant  Channel  Width  (cm) 


Radius  of  Fuel  region  (cm) 


8.532 

0.01 

0.25 

0.192 

0.15 

0.1 

0.15 

1.0 

1.205 


MATERIALS 

Emitter  Material 
Collector  Material 
Insulation  Material 
Cladding  Matoial 
Coolant 
Fuel  material 


W 

Mo 

AI2O3 

St^nless-  Steel 

NaK(78%) 

UO2 


OPERATING  CONDITIONS 

Coolant  Inlet  Temperature  (K) 

Coolant  Mass  Flow  Rate  (g/s) 

Load  Resistance  (mO) 

For  Steady-State  Analysis 
For  Transient  Analysis 
Initial  Steady-State  Fission  Power  (W) 

Nomiiud  Cs  Reservoir  Temperature  (K) 

REACTIVITY  FEEDBACK  COEFFICIENTS 

Fuel  Temperature  Reactivity  Feedback  Coefficient 

Emitter  &  Collector  Temperature  Reactivity  Feedback  Coeffici«it 

Tnailator  Temperature  Reactivity  Feedback  Coefficient 

Cladding  Temperature  Reactivity  Feedback  CkiefBcient 

Coolant  Temperature  Reactivity  Feedback  Coefficient 

Fuel  Doppler  Reactivity  Feedb^k  Coefficient 


850 

55.0 


5.0 

2.5 

1000 

567 


-1.2x10-5 

-6.1x10-6 

-8.7x10-6 

-1.0x10-5 

-2.0x10-6 

-0.8x10-6 


FISSION  THERMAL  POWER  (W) 


Figure  4.  Steady-stale  Results  Showing  the  Effect  of  Changing  die  Css  Pressure  and  the  Gap  Size 
on  die  TFE  Performance. 


CONVERSION  EFFICENCY  (%) 


fiarfon  powers  larger  than  1 150  W.  These  results  are  also  true  for  the  larger  gap  size  (0.25  nun)  at 
fission  powers  in  excess  of  1150  W.  At  lower  fission  powers,  however,  increasing  the  Cs 
pressure  in  the  larger  gap  results  in  ^ghtly  lower  conversion  efficiency  and  lower  electric  power  to 
file  load. 

Results  in  Fig.  4  show  fiiat  the  performance  parameters  of  the  TFE  at  a  low  and  hi^  Cs 
pressures  are  almost  identical  at  fission  powers  below  1150  W,  but  are  significantly  different  at 
higher  fission  powers.  At  a  low  fissicm  power  (  <  1150  W),  increasing  the  Cs  pressure  negligibly 
affects  the  operation  of  the  TFE,  while  increasing  file  gt^  size  always  results  in  both  hi^er  electric 
power  to  file  load  and  hi^er  converaon  efficiMicy.  For  example,  at  a  Cs  pressure  of  4  torr  the 
load  electric  power  increases  mtxiotcuiically  with  fission  power,  and  increasing  the  gap  size  from 
0.094  mm  to  0.25  mm  increases  the  electric  power  to  fiie  load  by  about  15%  and  50%  at  a  fission 
power  of  1400  W  and  800  W,  respectively.  Conversely,  at  the  low  Cs  pressure  (1.5  torr),  both 
the  load  electric  power  and  ccmversion  efficiency  increase  with  fission  power  until  fiiey  reach 
maTimiim  values,  then  fiiey  decrease  as  file  fission  power  is  increased.  Note  fiiat  decreasing  fiie 
gap  giy/"  decreases  fiie  maTimam  conversion  efficiency  and  the  maximum  load  electric  power,  but 
increases  the  conespcmding  fission  powers.  For  example,  at  a  Cs  pressure  of  1.5  torr  and  gap  aze 
of  0.25mm,  bofii  the  maximum  efficiency  ( - 10.8%)  and  the  maximum  electric  load  ( ~  122  W) 
occur  at  a  fissim  power  of  1150  W.  However,  i»1ien  the  gap  size  is  reduced  to  0.094  mm,  file 
maximum  conversion  efficiency  and  load  electric  power  decrease  to  9.2%  and  106  W, 
respectively,  but  file  corresptuiding  fisaon  power  increases  to  1200  W. 

hi  summary  ,  the  results  in  Fig.  4  suggea  fiiat  increasing  the  gap  size  always  enhances  the 
performance  of  the  TFE,  vtiiile  increasing  file  Cs  pressure  in  the  interelectrode  gap  inagnificantly 
affects  file  TFE  performance  at  low  fission  powers,  but  markedly  oihance  its  performance  at  high 
fission  powers.  However,  for  TFE  load-following  consideration  it  is  desirable  to  operate  at 
fissions  powers  below  that  corre^nding  to  fiie  maximum  load  electric  power.  As  such  fission 
powers,  it  is  advisable  to  lower  the  Cs  pressure,  hence  reducing  its  consumptitm,  with  little  effect 
cm  the  TFE  performance.  On  fiie  ofiier  hand,  the  results  clearly  demonstrate  that  a  reduction  in  the 
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gap  to  swelling  of  nuclear  fuel  would  strongly  degrade  the  performance  of  the  TFE. 

Hgures  5  and  6  present  flie  results  on  the  effect  of  fission  diermal  power  on  the  the  steady- 
state  temperatures  of  tiie  UO2  fuel,  W-cmitter,  and  the  Mo-coUector.  These  Hgures  indicate  the 
steep  rise  in  bofli  the  fuel  and  emitter  temperatures  as  the  fissimi  power  of  die  TFE  is  increased. 
Fbr  example,  increasing  the  fission  power  from  800  W  to  1400  W,  causes  the  fuel  and  emitter 
temperature  to  increase  from  approximately  2200  k  to  3050  K,  and  from  1700  K  to  1920  K, 
respectively.  Also,  the  coire^Kxiding  collector  temperature  increases  from  945  K  to  1020  K  [see 
Fig.  6].  Results  in  Rgs.  5  and  6  show  that  increasing  either  the  gap  oze  or  the  Cs  pressure 
insignificantly  affects  the  temperatures  of  the  TFE. 

4.2  Effect  of  Coolant  Temperature 

The  results  on  the  effect  of  the  coolant  inlet  temperature  on  the  TFE  conversion  efficiency 
and  the  load  electric  power  are  shown  in  Fig.  7.  As  this  figure  demonstrates,  the  inlet  temperature 
of  the  NaK  coolant  strmgly  affects  both  the  overall  efficiency  of  the  TFE  and  the  electric  power  to 
die  load,  paiticularly  at  hi^  fisrion  power.  For  a  gap  size  of  0.25  mm,  Cs  pressure  of  1.5  torr, 
and  fissinn  powers  in  excess  of  1100  W,  decreasing  the  coolant  inlet  temperature  causes  both  the 
conversion  efficiency  and  load  electric  power  to  increase.  Conversely,  at  lower  fission  powers, 
decreasing  the  inlet  coolant  temperature  causes  both  the  coiversion  efficiency  and  load  electric 
power  to  decrease. 

fr  is  worth  noting  that  the  fission  power  corresponding  to  the  maximum  load  electric  power 
("  1200  W)  is  independent  of  the  coolant  inlet  temperature,  while  that  corresponding  to  die 
iwaTimnin  c(Miversi(Mi  efficiency  increases  as  the  coolant  inlet  temperature  is  decreased.  Figure  7 
diows  that  at  a  fisrion  power  of  1400  W,  lowering  the  coolant  inlet  temperature  to  die  TFE  from 
850  K  to  800  K,  causes  both  die  TFE  conversion  efficiency  and  the  electric  power  to  the  load  to 
increase  by  about  12%,  from  ~  8%  to  9%  and  from  112.5  W  to  126  W,  respectively.  However,  in 
die  load-following  region,  at  a  fission  power  of  800  W,  lowering  the  inlet  coolant  temperature 
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Hguie  6.  Steady-state  Results  Showing  the  Effect  of  Changing  iht  Cs  Pressure  and  Gq>  Size  on 
die  Collector  Temperature. 
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Hguze  7.  Steady-state  Results  Showing  the  Effect  of  (Slanging  the  Coolant  Temperature  on  the 
TFE  Peifonnance 
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from  850  K  to  800  K,  causes  both  the  TFE  conversion  efficiency  and  the  electric  power  to  the  load 
to  decrease  by  about  4.5  %,  from  "  9.2  %  to  8.8%  and  from  73  W  to  70  W,  respectively. 

In  summary,  the  coolant  inlet  temperature  could  be  used,  in  addition  to  the  Cs  pressure,  to 
regulate  the  operatitMi  of  the  TFEs.  When  tolerating  at  hi^  fission  power,  niiere  the  TFE  is  non¬ 
load  following,  lowering  the  inlet  coolant  temperature  and/or  increasing  the  Cs  pressure  would 
inqirove  the  TFE  performance.  Conversely,  wh«i  operating  at  low  fission  power,  whwe  die  TFE 
is  load-following,  die  coolant  inlet  temperature  may  be  increased  to  improve  the  performance  of  the 
TFE. 

4.3  Transient  Operation 

Figures  8(a)  and  8(b)  present  the  changes  in  the  total  reactivity,  fission  power,  and 
temperatures  of  the  fuel,  emitter  and  collector,  following  a  $0.5  positive  step  input  in  reacdvity.  bi 
these  Figures  the  initial  steady-state  fission  thermal  power  before  the  transient  is  1,000  W  [see  Fig. 
8(b)].  Rgure  8(b)  shows  the  fission  power  and  the  temperatures  of  the  TEF  increasing  initially 
very  rapidly  due  to  the  step  increase  in  reactivity  and  the  rapid  increase  in  fission  power.  The 
negative  reactivity  feed-back  due  to  the  increase  in  temperature  progressively  slows  the  rate  of 
fiaginn  power  increase.  Eventually,  the  fission  power  as  well  as  the  temperatures  in  the  different 
regions  of  the  TFE  peak,  as  the  negative  reactivity  equals  the  input  reactivity  [see  Fig.  8(a)]. 
Beyond  this  point,  the  total  reactivity  becomes  negative,  forcing  the  fission  power  and 
temperatures  of  the  TFE  to  decrease.  Ultimately,  steady-state  is  reached  as  the  net  reactivity  in  the 
nuclear  reactor  core  becomes  equal  to  zero,  ~  250  seconds  after  the  initiation  of  the  transient. 

Rgures  9(a)  and  9(b)  present  tire  results  of  a  negative  step  input  of  reactivity  of  $0.5.  A 
comparison  of  the  results  in  Rgs.  8(b)  and  9(b)  indicate  that  foUowing  a  positive  and  a  negative 
reactivity  input  of  $0.5  die  fission  power  changes  initially  by  as  much  as  +350%  and  -40%, 
restively.  This  similar,  but  opposite,  response  is  also  indicated  not  only  by  the  change  in  the 
emitter  temperature,  but  also  the  changes  in  the  load  electric  power  and  conversion  efficiency  of 
the  TFE  [see  Figs.  10(a)  and  10(b)].  Another  major  difference  in  the  response  of  the  TFE  to 
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Figure  8(a).  Total  Reactivity  Transient  Response  to  a  Positive  Step  Input  in  Reactivity  of  $  0.5 
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Figure  8  (b).  Transient  Response  of  Fission  Power  and  TFE  Temperatures  to  a  Positive  Step  Input 
in  Reactivity  of  $  0.5 
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Hgmt  9  (a).  Total  Reactivity  Transient  Response  to  a  Negative  Step  Input  in  Reactivity  of  $  0.5 
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Figure  9(b).  Transient  Respcmse  of  Fission  Power  and  TFE  Temperatures  to  a  Negative  Step  Input 
in  Reactivity  of  $  0.5 
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positive  anti  negstive  reactivity  inputs  of  equal  magnitude  is  that  the  second  peakmg  in  the  fission 
thermal  power  for  file  ioiiner  (1,150  W)  is  much  larger  fiian  fiiat  for  the  latter  (1,005  W). 

fii  file  case  of  a  positive  reactivity  input  of  $0.5,  because  the  initial  rise  in  fission  power  is 
accompanied  by  a  much  slower  rise  in  load  electric  power  (see  Fig.  10(a)],  the  conversion 
efficiency  of  the  TEF  initially  drops  precipitously  from  about  10%  to  a  minimum  of  3%  [see  Rg. 
10(b)].  The  tninimum  efficiency  corresprads  to  the  peak  fission  power  occurring  after  the 
initiation  of  the  tranaent  Beyond  this  point,  file  TFE  efficiency  begins  to  rise  as  the  fisatm  power 
decreases  and  file  load  electric  power  emtinues  to  increase,  reaching  a  higb^  steady-state  value  of 
approximately  11  percent;  a  net  increase  of  one  percent.  A  amilar,  but  opposite  re^mises  in  the 
load  electric  power  and  TFE  conversion  efficiency  are  diown  following  a  negative  reactivity  input 
of  $  0.5  [see  Rgs.  10(a)  and  10(b)]. 

It  is  woifii  noting  that  the  transient  lespcmse  of  the  TI  cmiversion  efficiency  is  different  from 
fhat  described  above  for  the  TFE  overall  conversion  efficiency;  the  former  is  based  on  the  thermal 
power  input  to  the  emitter,  \diile  the  latter  is  based  on  the  thermal  power  ^nerated  in  the  fuel 
region.  While  the  TFE  efficiency  initially  decreases  (or  increases)  following  a  positive  (or 
negative)  step  input  in  reactivity,  the  TI  efficiency  increases  (or  decreases)  dowly  with  time.  At 
steady-state  the  TI  efficiency  equals  that  of  file  TFE,  whereas  the  fission  power  becomes  equal  to 
the  thermal  power  input  to  the  emitter.  At  the  beginning  of  the  transient,  the  change  in  the  fissKHi 
power  mainly  changes  file  thermal  energy  storage  in  file  fuel,  and  hence  its  temperature.  The 
fiiermal  input  to  the  emitter  changes  very  dowly  with  time  because  of  the  low  thermal  diffiisivity  of 
the  UO2  fuel  and  the  presence  of  an  open  gap  between  the  fuel  pellets  and  the  emitter  [see  Rg.  1]. 

Rgure  1 1  presents  file  operating  cmiditioo  before  and  after  the  transients  on  a  plot  of  the  J-V 
of  the  TI  ctmverter.  B  can  be  seen  that  initially,  file  TI  converter  operates  in  file 
ignited  obstructed  mode.  As  shown  in  Rg.  11,  a  positive  step  ii^ut  in  reactivity  of  $0.5  moves 
file  operation  conditions  of  the  TI  converter  closo-  to  the  transition  point  ^ ),  resulting  in  a 
hitler  convertion  efficiency  and  higher  load  electric  power  [see  Rg.  10(b)  and  11].  Ccmversely,  a 
negative  step  input  in  reactivity  of  $0.5  moves  the  operation  conditims  of  the  TI  converter  away 
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Figure  10(b).  Compaiis(»  of  the  Transient  Response  of  TI  and  TFE  Conversion  Efficiencies  to  a 
Positive  and  Negative  Step  hiput  in  Reactivity  of  $  0.5 
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Hguit  11.  J-V  Qiaracteiistics  Blustrating  the  Operation  Modes  of  the  TFE  Following  a  Positive 
and  Negative  Step  hiput  in  Reactivi^  of  $  0.5 


from  the  transition  point  ),  hence  lowering  both  the  TI  conversion  efficiency  and  the  load 

#  electric  power. 

4.4  Load-Following  Characteristics 

^  For  an  initial  steady-state  fisaon  power  of  1,000  W,  Hg.  12  presents  the  load-following 

diaracteristics  of  die  TFE.  As  this  Figure  diows,  because  of  the  negative  temperature  reactivity 
coefficients  of  the  TFE  [see  Table  1],  the  nuclear  fission  power  is  always  load  following.  An 

^  increase  in  the  load  electric  demand  increases  electttm  cooling  of  the  emitter,  dierefore  lowering  its 

temperature  and  that  of  the  fiiel  region.  Such  decrease  in  temperatures  introduces  a  positive 
temperature  reactivity  causing  die  fission  power  to  increase  commoisurate  with  the  load  demand. 

#  Conversely,  a  decrease  in  load  demand  increases  die  fuel  and  emitter  temperatures,  hence  lowering 
fission  power  doe  to  the  negative  temperature  reactivity  feedback  effect. 

While  the  fission  power  is  always  load-following  the  integrated  TFE  is  only  partially  load- 

#  following.  As  demrmstrated  in  Fig.  12,  the  load  electric  power  for  the  base  case  TFE  increases 
with  increased  load  demand;  it  peaks  at  98  W  for  a  100%  load  demand.  A  further  increase  in  load 
demand  is  met  with  a  decline  in  die  load  electric  power.  Therefore,  the  TFE  being  analyzed 

#  (corresponding  load  resistance  is  2.5  mW)  is  load-following  below  100  %  load  demand  and  is  ntm 
load-following  at  higher  load  demands,  ft  diould  be  noted  that  changing  the  design  and  operating 
conditirms  of  the  TFE  would  shift  the  peak  load  electric  power  either  to  the  right  or  to  the  left  of  the 

#  point  of  100%  load  demand  in  Hg.  12. 

Figure  13  presents  an  operating  surface  of  the  load  following  characteristics  of  the  TFE. 
Because  of  die  reactivity  feedback  effects,  changes  in  load  demand  causes  both  the  load  electrical 

#  power  and  the  fission  power  to  change.  Hgure  13  diows  that  ftir  a  given  load  electric  power  diere 
are  multiple  combinations  of  fission  power  and  load  (femand.  The  contours  in  Fig.  13  indicate 
constant  load  electric  power  regions.  The  dark  region  in  Fig.  13  corresponds  to  the  maximum 

#  attainable  load  electric  powers  between  120  and  140  W  for  the  base  design  TFE.  The  results  in 
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Hguie  12.  Load-following  Characteristics  of  the  TFE  Showing  die  Effect  of  Changing  the  Load 
Demand  on  Bodi  die  Load  Electric  Power  and  Bssion  Powo* 
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Fig.  13  should  be  useful  to  the  reactor  operator  for  determining  the  suitable  combinations  of  load 
electric  power  and  fission  power  needed  to  re^Kuid  the  dianges  in  the  load  demand. 

Results  presented  in  Fig.  14  demraistrate  the  effect  of  dianging  the  Cs  vapor  pressure  in  die 
interelectrode  gap  and/or  die  width  of  the  gap  on  the  load-following  characteristics  and  die 
conversi(Hi  efficiency  of  the  TFE.  For  the  base  load  demand  and  for  a  gap  width  of  0.25  mm, 
increasing  the  Cs  vrqior  pressure  from  1.5  torr  to  4.0  torr  insignificandy  affects  the  load-following 
response  and  the  conversion  efficiency  of  the  TFE,  since  it  does  not  change  the  operating 
condidcm  of  the  Tl  converter  [see  the  J-V  characteristics  in  Fig.  11].  However,  for  a  smaller  gap 
of  0.094  mm,  increasing  the  Cs  vapor  pressure  to  4.0  torr  diifts  the  operating  c(Hiditi(ms  of  the  TI 
converter  closer  to  the  )  transition  point,  thus  increasing  both  the  electric  current  and  the 
output  voltage,  resulting  in  a  higher  load  electric  power  and  a  slightly  higher  efficiency.  It  diould 
be  noted  that  at  high  Cs  pressure  the  power  losses  in  the  electric  leads  increases  due  to  the  high 
current  density,  causing  the  the  load  voltage  to  become  very  small  [see  Fig.  15].  Fig.  16(a)  shows 
die  effect  of  changing  the  Cs  vapor  pressure  and/or  the  width  of  the  interelectrode  gap  (« the  TFE 
fission  power  and  the  emitter  temperature,  while  Fig.  16(b)  presents  the  corresponding  effects 
die  load  current  and  voltage. 

These  results  suggest  that  in  a  TI  nuclear  power  system  it  is  desirable  at  the  beginning-of-life 
to  conserve  Cs  by  lowering  its  vapor  pressure,  since  increasing  the  Cs  pressure  would  not  affect 
die  load  electric  load  of  the  TFEs.  However,  diould  fuel  swelling,  after  operating  for  an  extended 
period  of  time,  decreases  die  width  of  the  interelectrode  gap,  the  load  electric  power  will  decrease, 
fii  this  case,  the  load  electric  power  could  be  restored  increasing  the  fission  power,  and  tally 
partially  by  increasing  die  Cs  vapor  pressure  in  the  gap  and/or  the  coolam  inlet  temperature. 

5.  SUMMARY  AND  CONCLUSIONS 

A  dynamic  model  that  simulates  both  transient  and  steady-state  operations  of  fully  integrated 
single-cell  TFEs  is  developed  and  used  to  perform  parametric  analyses.  The  steady-state  analytis 
investigated  the  effects  of  changing  the  Cs  vt^r  pressure,  tte  width  of  the  interelectrode  gap,  and 

38 


CONVERSION  EFRCIENCY 


LOAD  CURRENT  DENSITY  (A/cm*) 


Hgure  15.  Effect  of  Cs  ftessure  and  Gap  Size  on  the  J*  V  Characteristics  and  Operation  of  the 
TFE. 
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Hguie  16(a).  Effect  of  Cs  Pressure  and  Gap  Size  on  the  Emitter  Temperature  and  Fissimi  Power 
as  FUr  titms  of  Load  Demand. 
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Figure  16(b).  Effect  of  Cs  Pressure  and  Gap  Size  on  the  Load  Voltage  and  Churent  as  Fbnctions  of 
Load  Demand. 
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die  coolant  temperature  on  the  performance  as  well  as  the  load-following  characteristics  of  the 
TFE.  bi  addititm  to  diese  effects,  the  transient  analysis  investigated  the  response  of  the  TFE  to  a 
step  input  in  reactivity.  The  calculated  parameters  include  fission  power,  load  electric  power,  TFE 
conversion  efficiency,  h>^  current  and  volu^,  and  temperatures  in  the  different  regions  of  die 
TFE  (fuel,  emitter,  collector,  electric  insulator,  cladding  and  coolant)  as  functions  of  load  demand. 

Steady-state  results  diow  that  the  coolant  inlet  temperature  mig^t  be  used,  in  addition  to  the 
Cs  pressure,  to  regulate  the  operation  of  die  TFEs.  When  operating  at  a  high  fissitm  power,  udiere 
die  TFE  becmnes  ntm-load  following,  it  is  advisable  to  lower  the  inlet  coolant  temperature  in  order 
to  improve  the  TFE  performance.  Conversely,  vdien  operating  at  low  fission  power,  where  the 
TFE  is  load-following,  increasing  the  coolant  inlet  temperature  enhances  die  performance  of  die 
TFE.  b  actual  TI  power  system  appUcatims,  the  inlet  coolant  temperature  to  the  TFEs  could  be 
regulated  by  controllmg  the  the  mput  power  to  the  electromagnedc  pumps  vdiich  circulate  the  liquid 
metal  coolant 

Steady-state  results  also  show  that  mcreasing  the  gap  size  always  enhances  the  performance 
of  die  TFE,  while  bcreasmg  the  Cs  pressure  b  the  mterelectrode  gtq>  insignificandy  affects  die 
TFE  performance  at  low  fission  powers,  but  markedly  enhance  its  performance  at  bgh  fission 
powers.  However,  smce  its  desirable  to  operate  the  TFE  is  the  load-followmg  regime,  it  is 
advisable  to  lower  the  Cs  pressure,  and  hence  ledudng  its  consumpticm,  widi  little  effect  on  the 
TFE  operadcm.  On  d^  other  hand,  the  results  clearly  demonstrate  that  a  reducdm  b  the  gap  size 
due  to  sweUbg  of  nuclear  fuel  would  strongly  degrade  the  performance  of  the  TFE. 

Transient  results  show  diat  although  a  nuclear  reactor  widi  negative  temperature  reactivity 
coefficients  is  always  load  foUowbg,  die  TFEs  are  only  partially  load-foUowbg.  For  the  TFE 
base  case  design  and  operation  parameter,  the  load  electric  power  peaks  at  1(X)%  of  load  demand, 
which  coire^xmds  to  a  load  reastance  of  2.5  mW.  Beyond  this  pobt,  bcreasbg  die  load  demand 
towers  the  TFE  electric  power  power  ouqiut  (non  load-followbg  regime),  while  reducbg  the  load 
demand  is  acconqianied  by  a  reduction  b  the  load  electric  power  (load-followbg  regime).  Results 
also  diow  that  for  TFEs  with  large  bterelectrode  g^s,  it  is  desirable  at  the  beginnbg-of-life  to 
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conserve  Cs  by  lowering  its  vapor  pressure;  increasing  the  Cs  pressure  would  not  affect  the  load 
electric  load.  However,  should  fuel  swelling,  after  operating  for  an  extended  period  of  time, 
causes  the  width  of  die  interelectrode  g^  to  decrease,  both  the  TFE  ctmveraon  efficiency  and  the 
load  electric  power  will  decrease.  In  diis  case,  the  load  electric  power  mig^t  be  restored 
increasing  the  fission  power,  and  only  partially  by  increasing  the  Cs  vapor  pressure  in  the 
interelectrode  and/or  regulating  die  coolant  inlet  temperature. 
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